#########################################################################
###########################CARREGAR PACOTES##############################
#########################################################################
require(RcmdrMisc) #pacote para carregar arquivos do Excel
require(agricolae)
require(dae)


#########################################################################
#######################Planejamento do DQL###############################
#########################################################################
t = 5 ##numero de tratamento
trat=seq(1,t,1)
DQL=design.lsd(trat,  serie = 1, seed = 12345)
DQL1=DQL$book

Croqui1=matrix(DQL1[,4],t,t)
Croqui1

#########################################################################
############################Análise do DQL###############################
#########################################################################

#Mudar diretório
setwd("D:/OneDrive/UFMT/Disciplinas/Veterinaria/experimental/2019/dados")
setwd("C:/Users/ander/OneDrive/UFMT/Disciplinas/Veterinaria/experimental/2019/dados")

source("funcoes.txt")


##Importar dados
dados=readXL("anestesico.xls",sheet = 1,stringsAsFactors=TRUE)
attach(dados)
str(dados)

attach(dados)
str(dados)
linha=as.factor(linha)
tapply(trat,list(linha,coluna),sort)

######Análise de Variância#####################################3
model=aov(vr~coluna+linha+trat)
anova(model)
cv.model(model)

#####PRESUPOSIÇÕES DA ANOVA
#Gerar os erros padronizados do modelo
e=rstandard(model)



#Teste de Shapiro-Wilk - NORMALIDADE
shapiro.test(e) 

#Teste de Oneill-Mathews - Homogeneidade de Variância
VR=vr
Linha=linha
Coluna=coluna
Trat=trat
oneillmathewsDQL(VR,Trat,Linha,Coluna)

##Teste de Durbin-Watson - Independencia
set.seed(12345)
durbinWatsonTest(model)


###Teste de aditividade dos efeitos de tratamento, linha e coluna
model1=aov(VR~Coluna+Linha+Trat + Error(Linha*Coluna))
tukey.1df(model1, dados, error.term= "Linha:Coluna")